Fermi National Accelerator Laboratory 


>*7 

FERMILAB-Pub-9 1/2 4 0 -A 
September 1991 


Cosmological Constraints on Pseudo-Nambu-Goldstone Bosons 


Joshua A. FYieman 

NASA/Fermilab Astrophysics Center 
Fermi National Accelerator Laboratory 
Batavia , IL 60510-0500 

and 


Andrew H. Jaffe 

Department of Astronomy and Astrophysics 
The University of Chicago 
56^0 South Ellis Avenue, Chicago, IL 606S7 


Abstract 

Particle physics models with pseudo-Nambu-Goldstone bosons (PNGBs) axe characterized 
by two mass scales: a global spontaneous symmetry breaking scale / and a soft (explicit) 
symmetry breaking scale A. We investigate general model-insensitive constraints on this 
two-dimensional parameter space arising from the cosmological and astrophysical effects 
of PNGBs. In particular, we study constraints arising from vacuum misalignment and 
thermal production of PNGBs, topological defects, and the cosmological effects of PNGB 
decay products, as well as astrophysical constraints from stellar PNGB emission. Bounds 
on the Peccei-Quinn axion scale, 10 10 GeV < fpQ < 10 10 — 10 12 GeV , emerge as a special 
case, where the soft breaking scale is fixed at Aqcd — 100 MeV. 
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I. Introduction 


Particle physics models with spontaneously broken global symmetries are quite com- 
mon; by Goldstone’s theorem, the spectrum of such theories must contain a massless 
spin-0 boson for each broken symmetry group generator. If the global symmetry is only 
approximate, he., it is explicitly (in addition to being spontaneously) broken, the associ- 
ated bosons become massive pseudo-Nambu-Goldstone bosons (PNGBs). In nature, the 
best known example of this phenomenon is the 7r meson, associated with chiral symmetry 
breaking. Particle theory provides a host of additional PNGB candidates, 1 including ax- 
ions, majorons, familons, and schizons (although the majoron and familon may be exactly 
massless Nambu-Goldstone bosons). 

Such models are generally characterized by two mass scales: the spontaneous symmetry 
breaking scale /, and an explicit breaking scale A. The simplest example is that of a 
complex scalar $ with a potential V($*$) which is minimized at $ = /e**^//. The non-zero 
vacuum expectation value of $ spontaneously breaks the global 1/(1) symmetry $ — > $e m 
at the scale /, and the angular field <j) is the massless Nambu-Goldstone mode around the 
bottom of the ‘Mexican hat’ potential. At the lower scale A, a periodic potential for <f> 
of height ~ A 4 is generated. (The form of the potential may reflect a residual discrete 
symmetry.) The resulting PNGB has a mass given by m $ ~ A 2 //. In models with a 
large hierarchy between the scales / and A (/ ^> A), PNGBs are thus very light and 
also very weakly interacting, since their couplings are suppressed by inverse powers of /. 
Nevertheless, they can play an important role in astrophysics and cosmology. For example, 
in axion models, where / is the Peccei-Quinn scale /pq and the symmetry is explicitly 
broken by QCD instantons (through the chiral anomaly) at the chiral symmetry breaking 
scale Aqcd — 100 MeV, cosmological 2 and astrophysical 3 arguments constrain /pq to lie 
in a narrow window around /pq ~ 10 10 GeV (perhaps extending up to 10 12 GeV, although 
this point is controversial; see below). 

Although motivated by the strong CP problem, the QCD axion is a particular instance 
of a more general phenomenon. For example, Hill and Ross 4 have explored schizon models, 
in which A is associated with the mass of a light fermion (quark or lepton), and might 
plausibly lie in the eV (for neutrinos) to MeV (for charged leptons or quarks) range. 
On the other hand, superstring models contain one or more very heavy axion fields; for 
example, for the model-independent axion 5 (the imaginary counterpart to the dilaton), A 
is associated with the scale at which the gauge coupling in the hidden sector group becomes 
large. This naturally happens at a very high energy scale, 6 A ~ 10 14 — 10 17 GeV. (Indeed, 
it has been suggested that a PNGB with / ~ m p ] and A ~ 10 15 GeV is a natural candidate 
for the inflaton field. 7 ) 
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The lesson we draw from this is that there are a number of models with axions or 
axion-like particles in which the scale A has no connection with the QCD scale; indeed, 
from above we see that it may vary from the atomic to the Planck scale. We are thus led 
to consider the more general phenomenon of PNGB models in which both / and A axe 
a priori unconstrained. In this paper, we investigate the cosmological and astrophysical 
constraints on this two-dimensional parameter space. These bounds have previously been 
studied in detail for the QCD axion with fixed A = Aqcd; here, we explore how these 
constraints axe altered when the scale A is allowed to vary over a wide range. In Sec. II, we 
discuss PNGB production from initial vacuum angle misalignment and the bound which 
results from requiring that the PNGB density satisfy < 1. In Sec. Ill, we consider 

PNGB production from topological defects (strings, textures, and global monopoles) and, 
more generally, from long wavelength spatial gradients arising from the finite correlation 
length of the scalar field. We consider thermal production of PNGBs in Sec. IV. In Sec. V, 
we study a new wrinkle in PNGB phenomenology not present in ‘invisible’ axion models: 
if they axe sufficiently massive, PNGBs decay on a timescale shorter than the age of the 
Universe; we examine the concomitant constraints on their decay products. In Sec. VI, we 
discuss astrophysical bounds arising from PNGB emission from red giants and supernova 
1987a, and we conclude in Sec. VII. In all cases, the bounds we obtain axe relatively model- 
insensitive: aside from order unity couplings, the form of the PNGB interactions are fixed 
by symmetry, so they are not strongly model-dependent. 

Here, we briefly define some notation. In general, the PNGB <f> will be taken to have 
a (low-temperature) mass given by 


A 2 

= — = 10 5 eV 



( 1 . 1 ) 


Eqn.(l.l) may be thought of as defining the scale A in terms of the PNGB mass and /. 
In general, the PNGB field is related to an angular parameter of the broken symmetry, 
< i> = fO, where 6 runs from 0 to 2ir. we have absorbed any “winding” of the field into /. 


II. Non-Thermal Misalignment Production 

At high temperatures (T A), the angular degree of freedom representing the PNGB 
is randomly oriented. As the universe cools below a critical temperature T \ , the <f> field 
will roll to the minimum of its potential and will begin to oscillate coherently, resulting in 
a Bose condensate that can be treated as a classical field configuration. 

Expanded about the minimum of the potential, the Lagrangian density of the PNGB 
can be written as 

c = \<t> 2 ~ ( 2 . 1 ) 
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for a spatially homogeneous field <f>, to quadratic order in the field. (We consider the effects 
of inhomogeneities in the field in the next section.) In general, the potential will have 
higher order terms, but these are comparatively small for <j> < f (8 < 1); these anharmonic 
effects in the potential will slightly increase the PNGB energy density 8 if 6 > 1. 

We will parametrize the temperature-dependent mass of a general PNGB as 

m*(T) = amt(T = 0) (^) T > A (2.2) 

where a and v are fitting parameters of order unity. In the case of the QCD axion, 
the PNGB gains a mass at finite temperature due to QCD instantons 9 ; in this case, the 
behavior is approximately 8 ot ~ 0.1 and v 3.7. However, the form of Eq. (2.2) holds 
more generally; the temperature-dependent mass in schizon models has also been discussed 
recently 10 . 

The classical equation of motion for <f> is as usual for a scalar field in an expanding 
spacetime (ignoring for now the decay width of 

4> + 3Hj> + m^T) 2 ^ = 0. (2.3) 

At high temperatures where m^(T) <C 3 if(T), <f> = <f > i = fd\ is a constant and frozen to 
its initial value. As the temperature decreases, the mass turns on adiabatically; the field 
begins to oscillate coherently at a temperature T\ given by m^(Ti) = 3 H{T\). The PNGB 
energy density, averaged over an oscillation, then scales as 

H = ‘m.cn/tf, (2.4) 

where R(t) is the scale factor of the universe; thus the number of PNGBs per comoving 
volume is conserved, = p^/m^T) oc R~ 3 . [At low temperatures, when m^T) = 
m^( 0) = const., the PNGB energy density scales like non-relativistic matter, oc R 3 
(T < A).] As long as the expansion of the universe is adiabatic, the entropy density s is 
also proportional to R~ 3 , so the ratio n^/s has been constant since the onset of coherent 
PNGB oscillations. If the entropy in a comoving volume has increased by a factor 7 since 
the onset of PNGB oscillations, the ratio n^/s will simply be reduced by this factor. Thus 
we can find the present PNGB abundance by calculating the PNGB-to-entropy ratio when 
coherent oscillations begin; since p(2i) = m^(T 1 )^> 2 /2, we have 

n <i> ^ (2.5) 

s t x ~ 27t 2 0*iT 1 3 /45 

where s = 27r 2 p*T 3 /45 is the entropy density, g*(T) counts the number of relativistic 
degrees of freedom at temperature T, and the subscript 1 denotes the value at Ti . 
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Reexpressing this in terms of the initial misalignment angle 6 \ and the low-temperature 
mass m^, using ro^T)) = ZH\ = 5gl( 2 T? /m p i, appropriate for the radiation-dominated 
early universe, we can calculate the present PNGB density in terms of the critical density 
(defined by Hi = (8 ttG /3)p cr u) 

(n*/s)(s 0 /7W 

Pcrit Pent a \ 

a402 -1 ( 2 * 6 ) 

= 1.3 x 10- 10 h" 2 GeV " 2 


Here, the present Hubble parameter is H 0 = 100/ikm/sec/Mpc, and observations indicate 
0.4 < h < 1; we have also used the fact that the present entropy density is so = 2970 cm -3 
(for a photon temperature 11 of 2.735 K). Finally, T\ is given by ra^(Ti) = ZH(T \ ); for 
reference, in a radiation-dominated universe, 


Ti 

A 


Sjtff ) 

am v \ J 



Substituting Eq. (2.7) into (2.6), we find 


= 1.3 x 1O- 1O 0? 



7 -‘ K .?' +4 


(2.7) 


( 2 . 8 ) 


This expression merits several comments. First, it retains some implicit dependence on T\ 
through the term in g *\ ; however, between T ~ 1 TeV and today, g* has only changed by 
about two orders of magnitude, so this additional parameter does not have a substantial 
effect. Second, we have assumed in Eqns. (2.6) and (2.7) that the universe is radiation- 
dominated at T\. Thus, Eqn. (2.8) is only strictly valid if > 3 H(T eq ) ~ 10 -26 eV, 
where T eq is the temperature when the universe first becomes matter-dominated. Third, we 
note that for the axion, the most well-studied and well-understood PNGB, there is a theo- 
retical uncertainty in this value of about 8 10 ±0 ' 4 , coming mostly from the parametrization 
of the finite temperature mass; for the general case we might expect this to be larger. 

Solving this for A, and assuming from here on negligible entropy production, we find 


A 

GeV 


(2.3 x 10 9 )(2.4 x lO 18 )^! 



(2.9) 

We can use our knowledge of the expansion and age of the universe to restrict these 
parameters: since the Hubble parameter h > 0.4 and the age of the universe to > 10 10 yr 
(from globular cluster and nucleocosmochronology dating) we have the constraint on the 
cosmic density, 12 flo h 2 < 1. This requirement translates into an allowed region of the A-/ 
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plane, for fixed values of the 0 ( 1 ) parameters (i.e., a, v, <7*1, 7, $i). We show this region 
in Fig. 1. 

We have chosen as a typical value of the initial misalignment angle 6\ = 7r/\/3 — the 
RMS value of an angle randomly distributed between — 7r and n. On average, we expect 
this value to hold because the present Hubble volume comprises a very large number of 
regions that were causally disconnected at t\ — f(Tj ), within each of which was arbitrary. 
However, if the universe inflated during or after the PNGB symmetry breaking (at T ~ /) 
without reheating to temperatures above the symmetry-breaking scale /, we have no way 
of knowing the initial angle in our observable neighbourhood. In an inflationary scenario, 
it may be possible to escape this bound entirely if 6\ ~ 0 in our neighbourhood of the 
universe, since the observable Hubble volume would come from a single causally connected 
region within which 6\ was homogoneous. See refs. 13 for discussion of these issues. 

Figure 1 shows that this process alone excludes a large portion of the A-/ plane above 
the line representing Slh 2 = 1. Note the spread of 1-2 orders of magnitude by varying the 
model-dependent parameters about the “axionic” case. In particular, for small values of 
$i, the amount of excluded parameter space is decreased rapidly since is c< $1 • 

We have also shown, for reference, lines of constant mass m For the axion, located 
on the horizontal line at A = Aqcd — 100 MeV , misalignment production requires /pq < 
10 13 GeV , in agreement with previous analyses 2,8 . 

III. PNGB Production from Symmetry Breaking 

When the universe cools below T ~ / and the symmetry is spontaneously broken, 
global defects appropriate to the topology of the vacuum manifold will form 14 . In the 
simplest case (which includes axions), the relevent broken symmetry is 1/(1), the vacuum 
manifold S 1 has non-trivial first homotopy group, 7r 1 (5 1 ) = Z (the integers), and the 
defects produced are global strings; we first consider this case. 

At T ~ /, the massive complex field $ = rolls down to its minimum, |$j = 

a = j£, anywhere around the “mexican hat,” since the potential for the phase of $ (the 
PNGB field 4>) does not have discrete minima until T ~ A. Because the angular field <j> is 
initially uncorrelated on scales larger than the horizon, there will be closed loops in space 
around which <f> winds by 2 ttN y where N is an integer; by continuity, there must be some 
point within each surface spanned by such a loop where the angular field is undefined. 
These points are the cores of global strings, where the complex field is trapped in the 
false vacuum, $ = 0 . Due to the gradient and potential energy of the configuration, these 
strings have an energy per unit length fj, ~ 7 r / 2 In (fd), where d is a characteristic distance 
between strings. (The logarithm comes from the gradient energy in the field.) 

In the absence of inflation, the global string network should, as in the case of gauge 
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strings, quickly reach the so-called scaling solution 12 ; for long strings, the energy density 
then is of order 

^string ~ fit ~ 2 ^ Q ^ (3.1) 

Ptotal 

For the scaling solution to be maintained, almost all the energy in strings must be radiated 
away each Hubble time as PNGBs (this is the crucial difference between global strings and 
gauge strings, which primarily decay through gravitational radiation). Thus the change in 
the relative PNGB abundance over one Hubble time will be 

A(n*/s)~^- (3.2) 

where u is the average energy per radiated PNGB. This process will be effective from the 
epoch of spontaneous symmetry breaking (T ~ /) until the time when the PNGB mass 
m $ becomes comparable to the expansion rate H (at temperature T\ as before): the latter 
epoch corresponds to the time when: (1) the strings become connected by domain walls 
(at T ~ A ~ T \ ) and rapidly chop themselves up, and (2) the strings can no longer radiate 
into very low-frequency PNGB modes (since m $ ~ H(T\) ~ o>(fi)) which, at least in the 
Davis- Shellard scenario (see below), dominate PNGB production. Using the value of the 
Hubble parameter appropriate to a radiation dominated universe ( H ~ ~ T 2 /m pl ) 

and integrating Eqn.(3.2) gives 



For global strings, there is considerable controversy surrounding the value of a>(<), the 
average energy per radiated PNGB. Davis and Shellard 15 (DS) claim that the PNGBs 
should be radiated predominantly into low-frequency modes — wavelengths of order the 
horizon size (w ~ f _1 ) — but Harari and Sikivie 16 (HS) argue for a 1/1 spectrum, giving 
u> ~ In Rather than make a decision regarding these contentious issues, we will 
parametrize our ignorance: the density of particles produced by the decay of global strings 
is given by 



where the value of 5 is a result of this debate over the PNGB spectrum. For the DS 
scenario, 5 ~ ln(/U) (note that the logarithm comes from jJL rsj ~ f 2 ln(/*0); 

for the HS scenario, S ~ 1, since the In (ft) in ui (taken out of the integral since it is 
much more slowly varying than the rest of the integrand) cancels with that from \i. Since 
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Ti ~ A <C /, the energy density of string-produced PNGBs is 

A 4 


a^h 2 = 2.3 x 10~ n GeV ~ 2 S 




= 2.3 x 10~ n S 


(J± 

\G eV 2 / V°- 2am pJ 


)(^-) 

/ \0.2am p i/ 


* 7 +2 l 

j.T" 


(3.5) 


Aside from the factor of S, the abundance of string-produced PNGBs is comparable to that 
produced by the misalignment mechanism (compare Eqn.(2.8)), and may complement this 
mechanism or supplant it. Thus the string-produced density — and the resulting bound on 
/ and A — for the HS case is comparable to that for misalignment production. 

For the DS scenario, we need to compute the quantity In (/<i); using Eqn.(2.7) we find 
ft l = (3/ 2 /2aA 2 )(5 ff if 2 //am p i)-^ ( *' +2) , or 


!„(/<,) = 93 + 21„ 1„ (X) - In (^M’) (3.6) 


In solving the equation for A, we can neglect the In A term provided <C m p \, valid 
over most regimes of physical interest. Furthermore, for the axion with a ~ 0.1 and 
v ~ 3.7, the third term has a maximum value of about 2. Thus, we can approximate 
MM ~ 91 + [(i/ + 4)/(i/ + 2)] ln(//m p i) and we can solve for limits on A due to cosmic 
string decay: 


A 

GeV 


(5.1 x 10 10 )(2.4 x lO 18 )^ 


91 + (3f) W//10 12 GeV ) - 16] 


(3.7) 


for DS strings; remove the logarithmic factor in the denominator for the HS case. Despite 
the presence of the logarithm, the dependence of A upon / is still approximately that of a 
power law. If the DS analysis is correct, the PNGB abundance is increased by almost two 
orders of magnitude and the limits on the parameter space are correspondingly tighter. 

If the global symmetry group is not U( 1), other defects will form, but similar results will 
obtain as long as these defects also follow scaling solutions like Eqn.(3.1). For example, if 
the symmetry breaking is SU( 2) — ► 17(1), global monopoles will form. Each monopole has 
an energy ~ 47 T f 2 R, where R is a characteristic distance between monopoles. Numerical 
simulations 17 show that there should be 0(1 - 10) global monopoles per Hubble volume, 
R H~\ so the energy density in monopoles, p ~ (//t) 2 ? obeys the scaling solution. 

Global textures 18 which form, for example, in the symmetry breaking 0(4) — » 0(3) also 
obey such a solution, and the results above will obtain. In general, then, the defects arising 
in global spontaneous symmetry breaking give rise to a PNGB density comparable to that 
produced by the misalignment mechanism, up to factors of order ln(/fi) (which depend 
on details of the spectrum and field configuration for the defect). 
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The bulk of the energy density of global strings, monopoles, or textures is carried in 
the gradient energy of the PNGB field itself (unlike gauge strings or monopoles, in which 
the energy is carried in the core of the defect). Thus, even ignoring the defects themselves, 
“Kibble” gradients in the PNGB field (which arise from the field taking on different random 
values in initially causally disconnected regions) will also result in a population of PNGBs: 
as the horizon grows to encompass previously uncorrelated regions, the PNGB field will 
smooth itself, with a correlation length £ ~ 1 j H ~ t. At temperatures T > T \ , while the 
PNGB is still effectively massless, the gradient energy density is 

d ~ kv0 = 

1 mn (3 - 8) 

2 £2 J 

Once these long-wavelength fluctuations in <j> enter the horizon, we can consider them as 
a coherent state of <j> particles. Thus, spatial gradients in the PNGB field will produce a 
PNGB energy density again comparable to that produced by the misalignment mechanism. 

The excluded regions of the A-/ plane are shown in Fig. 2. In this figure and below we 
only show the line for one set of model- dependent parameters, those appropriate for the 
axion. Bear in mind that these lines will have a spread of 1-2 orders of magnitude (compare 
Fig. 1). We show results for both DS and HS strings; we expect that any PNGBs produced 
from either topological defects or spatial gradients should span this range. For the axion 
itself, at A = Aqcd — 100 MeV, HS strings restrict / < 10 12 GeV (not much different 
than the misalignment production bound), and DS strings restrict / < 4 x 10 10 GeV (an 
improvement of almost two orders of magnitude). 

When the soft symmetry breaking that gives the PNGB its mass occurs at T ~ A, 
further topological defects may form, and we must insure that they are cosmologically 
benign. In the ZV( 1 ) case, initially the PNGB field varies smoothly from 0 to 2i r around 
the string. When the symmetry is explicitly broken at T ~ A, the PNGB field is forced 
to its minimum, but the variation over 2 tt cannot be removed. Thus, domain walls will 
form, bounded by strings. If the PNGB potential has a unique (non-degenerate) minimum, 
each string is bounded by a single wall, and the string-wall system rapidly chops itself up 
and disappears. (This assumes the strings have not been inflated away.) If the PNGB 
potential has multiple minima, however, or if the strings are inflated away before the waills 
form, the domain walls would come to dominate the energy density of the universe and the 
microwave background would be strongly anisotropic 19 . Related phenomena may occur 
in other scenarios as well. For example, if the initial defects are global monopoles, the 
second symmetry breaking will result in monopoles connected by strings, leading to rapid 
monopole annihilation; in this case, the secondary defects are harmless. 
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IV. Thermal PNGB Production 

So far, we have assumed that PNGBs were not in thermal equilibrium in the early 
universe. That is, we have been assuming that microphysical scattering processes have 
had small impact on the abundance of PNGBs. Now, we consider the case in which 
PNGBs may couple to other constituents of the early universe and the conditions under 
which thermal production may be important. 

Since we wish to study constraints on PNGB properties that are as model-independent 
as possible, we do not know its couplings a priori. In general, however, the coupling 
strength will be suppressed by powers of (1//). We will consider two possible cases: 
(I) the PNGB is coupled to ordinary matter (e.g., quarks, leptons, and photons). (II) the 
PNGB is coupled only to the matter sector associated with the scale A, that is, to particles 
with masses m = (A/Aqcd ) m/(G), where m is the mass of a quark or hadron, and f(G) 
depends on the gauge group associated with the scale A. For example, a “techni-axion” 
coupled only to technicolor particles would fall in category (II), with A the scale at which 
the technicolor gauge group becomes strong (~ the electroweak scale). 

The interactions of the PNGB with photons and fermions are 1 

C = 9<t>-n<t>E • B + (4.1) 

fermions * 


with coupling constants 


a/27r 

94>~n = 9 7 

m f 

9<t>ff = 9f-j- - 


(4.2) 


Here, gg and g 1 are constants assumed to be of order one, although in some cases they 
may vanish ( e.g ., for the hadronic axion, which does not couple to electrons at tree level, 
g e = 0) or may be large (e.#., in Sikivie’s omion model 20 , g 7 ~ 10 5 ). In general, the values 
of these constants depend on the appropriate gauge-group charges, and to first order do 
not depend on / or A. Note that we have assumed the coupling to fermions is purely 
pseudoscalar, as is the case for axions, majorons, and familons in models without schizons. 

We will assume (following ref. 21) that the PNGB is produced in reactions like ab <— ► <pX 
in the early universe, with species a, b, and X all in thermal equilibrium. The PNGB will be 
produced in copious amounts by scattering processes if it was ever in thermal equilibrium 
t.e., if its rate of production, T, exceeds the expansion rate of the universe, H. In that 
case, it will have an abundance given by the equilibrium value for relativistic bosons, 
^EQ ^ 1 . 2 T 3 /tt 2 corresponding to a present energy density 


Q<t>,TEh 2 


10 

130 eV g*F 


(4.3) 
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where g„p is the value of g, at the PNGB freeze-out temperature Tp (defined by T(Tp) = 
H(Tp)). Here we have assumed the PNGBs are relativistic at freeze-out, Tp ra^; since 
we will see below that Tp ~ A, this just corresponds to the usual assumption that / ^> A. 

Thus the crucial model-dependence enters into the calculation of T/H to decide if the 
PNGB has been in thermal equilibrium. If the PNGB couples to normal matter (case 
(I) above), the calculation is the same as that for the axion, and only the value of / — not 
A — enters. The PNGB abundance is controlled by two processes. Before the quark-hadron 
transition at T ~ Aqcd ~ m/v/5 (where mjv — 1 GeV is the nucleon mass), the dominant 
production mechanism is PNGB photo- or gluon production in the presence of a heavy, but 
still relativistic, quark, Q (since the coupling strength is proportional to the fermion mass) 
Q'jf — ► Q<f> or QG — ► Q4>. After the quark-hadron transition, in the absence of free quarks, 
the dominant production mechanism is nucleon-pion scattering or pion-PNGB conversion, 
Nn -> N<f>. 

In the case of PNGB coupling only to particles with masses scaled by A (case (II) 
above), we will assume the same processes occuring for the “scaled” matter (and also 
a transition from scaled quarks to scaled hadrons at T ~ A): Qj — ► Q<f> and Nir — ► 
N<f>, where all fermion masses have been scaled by A ~ A/Aqcd> e -9 m Q = A rriQ ~ 
(A/Aqcd)toq- 

For PNGB photoproduction, the production rate T = tiq (<t|u|) ~ aT(mf / f) 2 , and 
H ~ T 2 /m p i in a radiation-dominated universe, so 



For pion-PNGB conversion, the nucleon abundance is Boltzmann-suppressed, tin — 
(mjvT) 3 / 2 e -m;v / r , and the cross-section (<r|v|) ~ T 2 / (fm n ) 2 , yielding 

^ ~ (m N T) 3 / 2 g 2 N T< A. (4.5) 

Due to the exponential decay of the nucleon abundance, the value of T/H will reach 
a maximum value very soon after the quark-hadron transition at T ~ A. When the 
PNGB couples to ordinary hadrons (case (I)), it is at that time in thermal equilibrium 

((I7#)t~a > 1) for 

/ < 6 x 10 *g N GeV (I). (4.6) 

We will define x — m^r/A; for normal matter couplings, x ~ 5. For scaled cou- 
plings (case (II)), x will differ; e.p., in an SU(N ) technicolor model, the technicolor 
scale A oc Aqcd / and the mass of the lightest technibaryon uiltb oc m^y/W so 
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x — (xqcd/Nqcd)-^ ~ 5iV/3. In terms of x, the rate for scaled matter is maximized when 

j ~ 5.2 x 10“^A (-^) x- 3 ' 2 e-‘. (4.7) 

where A = m/m = (A/Aqcd)(N/3) for SU(N) technicolor models. In this case, the 
requirement for thermal equilibrium becomes, defining y = x/5, 

/ < 6 x 10 8 y* GeV A 1 'V* /4 e“ 5(, '" 1)/2 (H) (4-8) 

or, for x = 5, N = 3, 

/ < 10 9 (A/ GeV yl 2 gf] GeV (II). (4.9) 

When these constraints are combined, the result is shown in Fig. 2. H^h 2 is greater 
than 1 for the region above the appropriate lines; the hatched marks denote the intersection 
of the regions for which the PNGB is in thermal equilibrium and for which fl^h 2 > 1. For 
the axion, which couples to normal matter, the constraint given by Eqn.(4.6) applies, and 
requires that / < 6 x 10 8 GeV for the axion to be in thermal in equilibrium. In order for 
thermal axions to make up a considerable fraction of the closure density, the value of / 
would have to be quite small, / < 5 x 10 4 GeV . Such a value of / would result in a very 
strongly coupled axion, already ruled out on astrophysical grounds (see section VI). 


V. Unstable Particles 

Thus far, we have assumed that the PNGB is stable (at least over times longer than the 
age of the universe). If it decays with a lifetime t into products which are still relativistic 
at the present time, the original PNGB density will have been redshifted away between 
the time of decay and today: 


R( r ) 

f^Decay Products (^o) ^ ^Stable 


where R is the scale factor of the universe, Rq is its present value, n is such that R oc t n at 
the time of decay, t ~ r, ftstabie is the density the species would have if it had not decayed, 
and flop is the present density of the decay particles. 

If the particle decays when the universe is matter-dominated (r > <eq = 4.4 x 
10 10 (f2 0 h 2 ) -2 sec ), n = 2/3, while if it decays in the radiation-dominated era (r < <eq)) 
n = 1/2. During the matter-dominated epoch, R(t)/Ro = 2.9 x 10 -12 (fioh 2 ) 1//3 (t/ sec) 2 / 3 . 
During the radiation-dominated era, R(t)/Ro = 2.4 x 10 -1 °y* 1 ^ 12 (t/ sec) 1 / 2 . Expressing 
this in more appropriate units, 


, R(r) _ f 2.0 x 1 CT 28 GeV 2 / 3 (f 2 0 /i 2 ) 1 / 3 r 2 / 3 r > *eq; 
U ' Ro 1.7 x 10 -22 G eV^V 1 / 12 T 1/2 r < t E Q ; 
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which we parametrize as 


_ a n 

m— — = a u t . 

it 0 


(5.3) 


On dimensional grounds the PNGB decay rate should go sis 1 T ~ m^// 2 , so the lifetime 
r = 1/r would be given by 

(5.4) 


(wbv 


y^y 


where gd is a dimensionless coefficient of order one. (For comparison, the present age of 
the universe to — /(fi)ifj" 1 = 4x 10 17 sec for = 1, h = 1/2.) This gives 


fi D p(<o) = A n gZf 5n A- 6n n s (to). 


(5.5) 


We can now insert our previous expressions for Sis, eqns. (2.8), (3.5), and (4.3), and 
these can then be solved for A. For example, PNGBs produced via the misalignment or 
defect mechanisms have a stable density given by equations (2.8) and (3.5): 

fy, mis = C v Af p ; p = (5-6) 

where C v is the coefficient of A f p in equations (2.8) or (3.5) and depends upon v, a, </*j , 
etc. Using this in our results for the density of decay products gives 


Slopf = .4 n C„ ! ,?/ 5 ” + > , A 1 " 6 ” 


(5.7) 


or, solving for A, 

a = {A; 1 c; 1 g j n r 5n - p nvph 2 ) r: ^ . ( 5 . 8 ) 

In order to convert this expression, and similar ones for the production of PNGBs via 
string decay and thermal mechanisms, into limits on the PNGB parameters, we need to 
find a limit on flop- Clearly, requiring Slvph? < 1 would be the most conservative bound. 
However, if we make some assumptions about the nature of these decay products and the 
growth of structure in the vmi verse, we can put far stricter bounds upon flop- 

Observations of the cosmic microwave background 11 have determined that its devia- 
tions from a perfect blackbody are remarkably small. Specifically, spectral distortions due 
to the cosmological Sunyaev-Zel’dovich mechanism are below the level of current measure- 
ments. When a batch of photons are dumped into the universe (due to a process like 
PNGB decay), the primary photons will Compton scatter off of free electrons, and the 
resulting hot electrons will in turn Compton scatter the microwave background, moving 
photons from the Rayleigh-Jeans to the Wien region. Clearly, this process is only relevant 
when the photons are produced before the era of recombination: (r < tree)- The amount 
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of distortion is quantified by the Compton y-parameter y = (l/4)(Ap 7 /p 7 ), where A p 7 is 
the excess photon density injected by the decays; COBE measurements 1 1 have found that 
y < 3 x 1CT 4 . (Actually, at temperatures T > 10 4 K, the spectral distortion takes the form 
of a non-zero chemical potential y, but the resulting bound is similar; we therefore focus 
on y.) Once the excess photons are dumped into the universe, they evolve the same way as 
the already-present background photons — y does not vary with time and is proportional 
to the ratio of flop to ft 7 . 

We parametrize the amount of photon energy dumped into the universe by PNGB 
decay as Ap 7 (r) = B y f c p^(r), with B~, the branching ratio for decay to photons, and f c the 
efficiency of Compton rescattering at the time of PNGB decay (/ c ~ 1 for T > T reC i fc — 0 
for T < T iec ). Using this for PNGBs with a lifetime r < t rec ~ 5.6 x 10 12 (fio^ 2 ) -1 ^ 2 860 ? 
we find 

(5 - 9) 

or, since f l y h 2 = 2.6 x 10 -5 , 

f^DP^ 2 ^ — ■ (5.10) 

"7 Jc 

Although this limit on the CMB distortion via the Sunyaev-Zel’dovich effect is well- 
understood, it is only effective (/ c ~ 1) for particles decaying at redshifts 10 7 > zd > z rec ~ 
10 3 *.e., with lifetimes 10 6 sec < r < t rec ; photons injected at earlier times, t < 10 6 sec , 
will thermalize and produce no distortions. This bound also relies upon assumptions about 
the PNGB decay branching ratio to photons, B 7 . Because of these constraints, this bound 
actually gives us no new information, since the area of parameter space with the proper 
PNGB lifetime does not overlap with the area that is newly allowed by eq. (5.10). 

However, we can get a bound on flop independent of the decay products, so long 
as they are relativistic, by considering the formation of structure 22 . We know that the 
universe is matter-dominated today (as long as there are no exotic relativistic particles 
which came to dominate at very recent epochs; this caveat includes the decay products of 
the PNGBs, so we cannot put this bound upon flop if t ~ to)? that it must have been 
matter-dominated for long enough to allow structure formation to occur, since only in a 
matter-dominated epoch do density perturbations grow sufficiently rapidly, linearly with 
the scale factor of the universe. The epoch of matter domination is given by 

1 + Z eq = Ro/ReQ = f^NR/flrel (5.11) 

where f2 re i is the density of relativistic particles, fl re i = ft 7 +„+f2DP) and ^nr is the present 
energy density of non-relativistic species. Thus, a significant flop will decrease the redshift 
of matter-radiation equality. Moreover, if Odp > fl 7 +t/, then the PNGBs dominated the 
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universe at the time of their decay, and the universe went through two periods of matter- 
domination: a first phase dominated by the PNGBs before they decayed into relativistic 
particles, and a second phase dominated by the presently non-relativistic matter. At the 
final epoch of matter-radiation equality, the comoving wavelength of the horizon scale is 
(assuming three light neutrino species) 

Kq = 5 Mpc (^NRh 2 ) -1 ^1 + ■ ( 512 ) 

If the non-relativistic matter today is cold dark matter, k ~ approximately sets the scale 
where the present fluctuation spectrum makes the transition from the primordial spectrum 
on large scales, |6* | 2 ~ k n for A: -1 >• k ~^ , to the processed spectrum on small scales, 1 2 ~ 

k n ~ 4 for k~ 1 < k~^ . Enlarging this scale by having an appreciable flop will increase the 
large scale power in 6p/p , for a fixed small-scale normalization (8p/ p « 1 at a scale of 
8 h~ l Mpc today, modulo biasing). The extra large scale power will in turn increase the 
quadrupole anisotropy of the CMB, which has been constrained by the Relikt experiment 23 
to be (AT/T ) rm3] ; = 2 < 1.5 x 10 -5 at the 95% confidence level. (This bound assumes 
a scale-invariant (Harrison-Zel’dovich) primordial perturbation spectrum; the spectrum- 
independent bound is roughly a factor of two higher.) 

Recently, the qualitatively similar effects of a decaying l7keV neutrino on structure 
formation 24 have been analyzed. For a primordial scale-invariant spectrum, the quadrupole 
bound limits its lifetime to be r < 10 yr; we can translate this constraint into the require- 
ment that 25 

Kq ~ 54 h -1 Mpc. (5.13) 

Using this in equation (5.12) with 17 nr ~ 1, we then obtain a limit upon flop: 

^Dph 2 < 115/i 2 fi 7+t ,/i 2 = 5 x 10- 3 fi 2 , (5.14) 

or, for h = 1/2, 

flDph 2 < 1.3 x 10 -3 . (5.15) 

This corresponds roughly to the requirement that the universe become matter domi- 
nated before or during the epoch of recombination. This constraint on flop from structure 
formation is shown in Fig. 3, for each of the production mechanisms examined above. 
Decay processes make new areas of the parameter space “cosmologically allowed”: as we 
move away from t = to towards shorter and shorter lifetimes (t.e., earlier and earlier de- 
cay times), the density of PNGB decay products decreases because it is redshifted away. 
Note that the region of allowed parameter space for thermally-produced decaying PNGBs, 
especially those that couple to normal matter, is extremely small; any change in the 0(1) 
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parameters that control their abundance is likely to make this region even smaller or 
eliminate it entirely. Furthermore, we have throughout assumed that each limit applies 
independently to each production mechanism. If, however, more than one mechanism is 
effective, then the limits will be tighter, due to separate contributions to SI. 

VI. Astrophysical Constraints 

As we have already mentioned, the strength of the coupling of the PNGB to matter and 
radiation does not depend on the energy scale A, but is merely proportional to 1//. In fact, 
any work that restricts the couplings of any pseudoscalar particle is directly applicable to 
our general PNGBs. Astrophysical methods have been extremely valuable in constraining 
the mass of one pseudoscalar, the axion, and we will generalize the most important of these 
results here. 

If a PNGB couples strongly enough to matter or radiation, then it will be produced 
in copious amounts in astrophysical objects such as ordinary stars, red giant stars, and 
supemovae-if it is also light enough that its production is not Boltzmann-suppressed ( m ^ 
less than a few times the temperature of the object). If its interactions are still sufficiently 
weak, however, it may be able to stream freely out of the star after its production. This 
“cooling” mechanism — a misnomer, since the necessity of maintaining an equilibrium con- 
figuration generally increases the temperature — may strongly affect the evolution of these 
objects. 

The strongest stellar-evolutionary bounds to pseudoscalar couplings come from the 
examination of helium-burning (horzontal branch) stars. The coupling to photons is con- 
strained by the lifetime of these stars, observed via the fraction of stars seen in that area of 
a star cluster’s color-magnitude diagram. In the presence of PNGB cooling via the PNGB 
Primakoff process, 7 + (Z, e“) — ► <j> + (Z, e _ ), the length of this evolutionary stage will be 
shortened; a value of 

g+y. r < 10 - 10 GeV~ l or / > lO 7 ^ GeV (6.1) 

is necessary to keep the helium burning lifetime within an order of magnitude of obser- 
vations. This constraint holds for PNGB masses <C T. If the PNGB mass is greater 
than the core temperature for typical horizontal branch stars T — 8.6 keV, the PNGB 
emission rate will be suppressed by roughly the factor exp (— m^/T); this ignores any other 
mass-dependence of the emission rates, which will be swamped by the Boltmann factor for 
rritf, T. (We have worked out the limits from Compton emission in horizontal branch 
stars using the exact temperature-dependence of the rate from Raffelt and Starkman 26 
and the constraints on / and A make only a small difference in fig. 4 below.) In general, 
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emission rates are proportional to #^ 77 , so the constraint will become 

94>-n< 10- 10 e m / 17 ' 2keV GeV- J or / > lO 7 e" m '/ 17 - 2keV 0 7 GeV (6.1a) 

valid for either m# <C 8.6 keV or 8.6 keV. 

The constraint on the coupling to electrons relies on the physics of helium ignition in 
low-mass red giants. These stars have cores of helium nuclei and degenerate electrons. 
Due to this degneracy the pressure does not increase as the temperature rises: when 
it rises sufficiently for helium burning (via the triple-alpha process) a thermal runaway 
occurs until the temperature is sufficient for degeneracy to be lifted. If PNGBs couple 
sufficiently strongly to the degenerate electron gas, they will be able to cool the core 
effeciently enough to prevent the ignition of helium — and this time the term “cooling” is 
appropriate as the temperature is decreased due to the PNGB streaming. The PNGBs axe 
produced largely via the process of PNGB bremsstrahlung, e~ + Ze — ► e~ -f Ze + <f> with 
an internal photon line. Calculations 26 show that this cooling will prevent helium burning 
f° r 9<t>ee 3 x 10 -13 . In order for the cooling to actually be effective, however, the PNGBs 
cannot couple so strongly that they remain trapped in the core, keeping the temperature 
high. This results in a “window” for which helium ignition will not occur: 

3 x 1CT 13 < g^ ee < 6 x 10 7 ; (6.2) 

that is, helium ignition will occur normally for 

/<10 3 <7 e GeV or />10 9 <7 e GeV. (6.3) 

In order for these bounds to hold, the mass of the PNGB must be much less than the 
temperature. For heavy PNGBs, the cooling limit is modified as in equation (6.1a). In 
the trapping regime, the primary effect of masses > T will be to slow the PNGBs 
to non-relativistic velocities, decreasing the effective coupling a# = < 7 ^ ee / 47 r by a factor of 
the average thermal velocity of the PNGBs, a ^ = cv(m^ ) /T') 1 , thus reducing the coupling 

constant g,f, ee by a factor of (m^/T) 1 / 4 . In order to interpolate between the limits of 
m <t> <C T and > T, we will modify this factor to (1 + m^/T) 1 / 4 . Thus our Helium 
ignition limits become 

/< 10 3 (1 + m <fr /8.6keV) 1 /VGeV or / > 10 9 e~ m */ 17 - 2keV ^ e GeV . (6.3a) 

The strictest astrophysics! bounds on pseudoscalar couplings arise from observations of 
the duration of the neutrino pulse from supernova 1987a. If the nucleon coupling is strong 
enough, PNGB emission will cool the supernova core so efficiently that the neutrino burst 
duration will be shortened. If they couple too strongly, however, PNGBs will be trapped 
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in a “PNGB sphere” (analagous to the neutrino-sphere) which will again lengthen the 
neutrino burst sufficiently to coincide with observation. The dominant emission process in 
the SN core is PNGB-Nucleon bremsstrahlung, NN -*■ NN + <f> with an internal pion line. 

Numerical and analytical studies 2 ** have shown that the neutrino pulse will be unac- 
ceptably shortened for 

10" 10 < g<t>NN < 2 x 10- 7 (6.4) 

assuming equal coupling to neutrons and protons (g^nn = g<t>pp = g<l>NN )• Thus the 
neutrino burst will occur normally for 

/ < 5 x 10 *g N GeV or / > 10 10 <7 N GeV . (6.5) 

If we take account of masses T co re ~ 30 MeV as above, these constraints become 

approximately 

/ < 5x 10 6 (l + m^/30MeV) 1/4 5ivGeV or /> 10 10 e _m * /60MeV gu GeV . (6.5a) 

Of course, none of these astrophysical bounds obtain unless the PNGB couples to 
normal matter and electromagnetic radiation. If, for example, the PNGB couples only to 
the photon, then only the bound on g * 77 from the helium burning lifetime will matter. 
Similarly, in the case of the hadronic axion, which does not couple at tree level to the 
electron (i.e., g t is very small) only the photon and nucleon bounds are relevent. 

We show the astrophysical bounds in Fig. 4; for modest PNGB masses, the limits are 
merely upon the parameter / that controls the PNGB couplings. For extremely heavy 
PNGBs, however, the masses are comparable with the temperatures of the astrophysical 
objects under consideration, and the limits are strongly curtailed. 

For small values of /, the coupling of the PNGB to matter becomes strong enough that 
it would show up in terrestrial accelerators. On the basis of nondetection of T + + 7 

and J/V» —></> + 7 events, we can limit 1 

/> 10 3 GeV. 

This bound holds for PNGB masses m $ < lMeV , above which the decay to electron- 
positron pairs becomes important. Thus, the lower left of fig. 4 is also disallowed, despite 
the regions allowed by trapping in helium-burning stars. 

Other groups have also been investigating various astrophysical constraints on low- 
mass PNGB properties. Gnedin and Meszaros 27 have examined the conversion of PNGBs 
to radio photons in the magnetic field of a pulsar. Engel et al. 28 have used the nondetection 
of a direct axion signal from SN87a by the Kamiokande II detector to exclude the range 
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9.4 xlO 2 GeV < f / 9N ~ 1.0 x 10 6 GeV , for <C 30MeV. Ressell 29 uses the nondetection 
of a spectral axion line in the night sky resulting from present day decay of clustered axions 
to exclude //j T < 2 x 10 6 GeV , for PNGBs with decay times r ~ to- 

VII. Conclusions 

We have considered the cosmological and astrophysical bounds on the parameter space 
of broken global symmetry models with pseudo-Nambu-Goldstone bosons. As Figs. 1 and 
2 indicate, as the soft breaking scale A is increased, the cosmological upper bound on 
the spontaneous symmetry breaking scale / (due to the requirement that PNGBs have 
less than closure density today) is generally reduced, in accordance with expectations. 
However, this general trend is subverted when A becomes large enough that the PNGB 
decays on a timescale significantly shorter than the age of the Universe. In that case, since 
the relativistic PNGB decay products redshift away, the bounds are somewhat relaxed, 
83 Fig. 3 shows. Since the stellar emission of PNGBs depends only on the scale /, the 
astrophysical lower limits on the axion scale f PQ (from red giants and SNl987a) can be for 
the most part taken over to the general case. However, these bounds do have an indirect 
dependence on the scale A through the Boltzmann suppression factor for PNGB emission 
~ since = A 2 //. Thus, for large values of A, the astrophysical lower bounds 

on f are also relaxed. From Fig. 5, which summarizes the astrophysical and cosmological 
constraints on the parameters A and /, we see that there are essentially two allowed regions 
of parameter space: (i) for A < 1 - 50 GeV (depending on the DS vs. HS debate) and 
/ ^ 10 10 GeV, there is a triangular region allowed for stable, light (m^ < 10“ 5 eV) PNGBs; 
here the upper bound on / due to the cosmic PNGB density increases as A is reduced, 
(ii) For A > 10 — 1000 GeV , there is a second allowed region; as A is increased above this 
limit, a widening window for / opens up. In this region, the PNGBs are heavy (m^ > 1 
MeV) and unstable on a timescale short compared to the age of the Universe. 
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FIGURES 


FIG. 1. The excluded area of the A-/ plane for non- thermal misalignment production 
of stable PNGBs. The region above these lines has ft 0 A 2 > 1 and is excluded. The 
line labeled “base” is the axionic case, with the following parameters: v = 3.7, a = 0.1, 
</*! = 100, 0i = 7 r/\/3. The other lines differ from this case as marked in the legend. In 
addition, lines of constant mass (m* = A 2 //) are shown for reference. 

FIG. 2. The excluded area of the A-/ plane for all production mechanisms of stable 
PNGBs. For Misalignment and String production, the area above the ft* ft 2 = 1 lines 
is excluded as in Fig. 1. For Thermal production, the excluded area is above the line 
marked “In TE” and to the left of either the line marked “Thermal (Normal)” or “Thermal 
(Scaled),” depending on the PNGB couplings (see text). In this figure and below, we have 
set g N - 1. 

FIG. 3. The excluded area of the A-/ plane for decaying PNGBs, limited by constraints 
of structure formation and CMB distortion, ftopA 2 ^ 1.3 X 10 3 . The lines are as in 
Fig. 2, with the constraints bending over as the density redshifts away for PNGBs that 
have decayed by the present time. In this figure, we have set g* = 1- 

FIG. 4. The excluded area of the A-/ plane due to astrophysical limits on PNGB cou- 
plings. The regions below the curves are disallowed by the various astrophysical arguments 
as marked. In this figure, we have set g e = gw — 9-y — 1 text). Lines of constant 
PNGB mass are shown for reference. 

FIG. 5. The allowed areas of the A-/ plane with all constraints from Figs. 3-4 
combined. Solid lines are constraints from misalignment and topological defects; dotted 
lines are from thermal production; short dashed lines are from astrophysical constraints. 
“PNGB dark matter” labels the region where ft*/i 2 ~ 1. Lines of constant PNGB mass 
are shown for reference (long dashes). 
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